Scattering of light by molecules of light 
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We study the scattering properties of optical dipole-mode vector solitons recently predicted theo- 
retically and generated in a laboratory. We demonstrate that such a radially asymmetric composite 
self-trapped state resembles "a molecule of light" which is extremely robust, survives a wide range 
of collisions, and displays new phenomena such as the transformation of a linear momentum into an 
angular momentum, etc. We present also experimental verifications of some of our predictions. 
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Since the begining of the history, physics has studied 
simple objects and the way they arrange to form more 
complex objects. Some remarkable success included the 
atomic theory of matter, the discovery of the structure 
of nucleus in terms of protons and neutrons and, more 
recently, the substructure of nucleons in terms of quarks. 
Such concepts seem to be restricted to "solid" objects 
only, and at a first sight it could seem that nothing sim- 
ilar is possible for light. However, elementary robust ob- 
jects made of light have been known since the 70 's. In 
fact, spatial optical solitons — self-trapped states of light 
with particle-like properties — have attracted a consid- 
erable attention during last years as possible building 
blocks of all-optical switching devices where light is used 
to guide and manipulate light itself JjJ. Recent progress 
in generating spatial optical solitons in various nonlinear 
bulk media allows to study truly two-dimensional self- 
trapping of light and different types of interactions of 
multi-dimensional solitary waves Q . 

Robust nature of spatial optical solitons as self-trapped 
states of light that they display in interactions , allows 
to draw an analogy with atomic physics treating spatial 
solitons as "atoms of light" . Furthermore, when several 
light beams generated by a coherent source are combined 
to produce a vector soliton, this process can be viewed as 
the formation of composite states or "molecules of light" . 

Recently, we have predicted theoretically the existence 
of a robust "molecule of light" , a dipole-mode vector soli- 
ton (or "dipole" , for simplicity) that originates from trap- 
ping of a dipole beam by an effective waveguide created 
by a mutually incoherent fundamental beam f| . The first 
observation of this novel type of optical vector soliton has 
been recently reported in Ref. B , where the dipoles have 
been generated using two different methods: a phase im- 
printing and a symmetry-breaking instability of a vortex- 
mode composite soliton [EJ. 

The concept of vector solitons as "molecules of light" 
should be compared with photonic microcavity struc- 
tures, micrometer-sized "photonic quantum dots" that 
confine photons in such a way that they act like electrons 
in an atom m. When two of these "photonic atoms" 



are linked together, they produce a "photonic molecule" 
whose optical modes bear a strong resemblance to the 
electronic states in a diatomic molecule like hydrogen [Q . 
Self-trapped states of light we study here can be viewed 
as the similar photonic structures where, however, the 
photonic trap and the beam it guides are both made of 
light creating self-trapped photonic atoms and molecules. 

In this Letter we study the scattering properties of 
dipole-mode vector solitons ("molecules of light") and 
analyze, in particular, the interaction between these ob- 
jects and other robust structures made of light: scalar 
solitons ("atoms of light") and other dipoles. We de- 
scribe a number of interesting effects observed in such 
interactions, e.g. the absorption of a soliton beam by 
a dipole and replacement of the soliton with a dipole 
component, transformation of a linear momentum into 
an angular momentum with subsequent dipole spiraling, 
etc. Additionally, we verify experimentally some of these 
predictions for composite spatial solitons generated in a 
photorefractive crystal. The versatility of phenomena de- 
scribed here makes the dipole-mode vector soliton a com- 
plex object with promising applications in integrated op- 
tics in addition to its fundamental interest. 

The model and dipole solitons. We consider the prop- 
agation of two coherent light beams interacting incoher- 
ently in a saturable nonlinear medium. In the paraxial 
approximation, the beam mutual interaction can be de- 
scribed by a system of two coupled nonlinear Schrodinger 
(NLS) equations §|§, 

du 1 . , . , „ 

i-Q- z = --A x u + F(I)u, (la) 

= -lA ± v + F(I)v, (lb) 
az 2 

where u(r±,z) and v(tj_,z) are dimensionless envelopes 
of the beams self-trapped in the cross-section plane rj_ = 
(x, y) and propagating along z. The function F{I) — 
7(1 + sl)^ 1 characterizes a saturable nonlinearity of the 
medium, where s is a dimensionless saturation parameter 
(0 < s < 1) and I — \u\ 2 + \v\ 2 is the total intensity. 

Equations (1) describe different types of spatially local- 
ized composite solutions. The dipole-mode vector soliton 
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(or "a molecule of light") is a stationary state which is 
composed of a nodeless beam in the v component and a 
dipole beam (or a pair of out-of-phase solitons) in the u 
component. Solitons in the u component have opposite 
phases and thus they tend to repel each other, but the 
role of the beam v is to stabilize the structure making it 
robust. A numerical analysis of the linearized equations 
(1) shows no signs of linear instability of this composite 
structure ||. Moreover, it was shown || that such robust 
dipole-mode vector solitons exist for a wide range of the 
beam powers P u = J \u\ 2 dr± and P v = j\v\ 2 r^. Since 
we are interested in showing stability far from the regime 
in which one beam is dominant and acts as a waveguide, 
all numerical experiments are performed using as initial 
states for the collisions experiments stationary states in 
which P u ~ P v . 

Qualitative analysis. We are interested in the dynam- 
ics of the dipole under the action of finite external per- 
turbations introduced by its collision with other objects. 
The word "finite" emphasizes the fact that we can no 
longer make use of linearized equations and that we must 
deal with the full system (1). This fact, combined with 
the complex structure of the dipole which lacks radial 
symmetry, makes analytical predictions on the dipole dy- 
namics very difficult. Nevertheless, as will be shown be- 
low, one may extract some general rules on which quali- 
tative predictions may be based. 

The idea is that the dipole can be seen as a bound 
state of a soliton beam (in u) plus a pair of vortices with 
opposite charges (in v), and therefore many effects ob- 
served in the beam collisions can be understood once the 
mutual interaction of these simpler objects is known. 

One of the components of the dipole is a soliton beam 
(to be referred to as soliton hereafter). Spatial solitons 
are stable localized states which have no nodes and which 
are the states of minimum energy of the system for a 
fixed power. When two of these solitons are in differ- 
ent beams (say, one in u and the other one, in v), they 
interact incoherently and attract each other. Thus, dur- 
ing an incoherent interaction two solitons attract each 
other and either become bound or scatter. In the former 
case, we have an example of what we call a molecule of 
light, which is typically referred to as "vector soliton" . 
However, when two solitons are in the same component, 
their mutual interaction depends on their phase differ- 
ence. When this quantity is small or zero, they interact 
attractively, whereas if their mutual phases differ by tt, 
they repel each other. 

Another nonlinear structure that should be mentioned 
in this context is a vortex-mode composite soliton intro- 
duced in H which in our model (1) is an unstable object 
(see details in ||). Vortices may only be partially stabi- 
lized (i.e., their decay rate becomes smaller) by sharing 
space with a large soliton beam (e.g., when a vortex in 
the beam u is guided by an effective waveguide created 
in the component v). Thus, a composite state of a vortex 



plus soliton constitutes an unstable molecule of light. 

Concerning a dipole, it can be seen as a pair of vortices 
as described above or, alternatively, as a bound state 
of two solitons with a phase difference of tt. Although, 
in principle, these solitons should repel each other, the 
system is stabilized due to the interaction with a soliton- 
induccd waveguide created in the other component. 

Soliton- dipole scattering. The first type of numerical 
simulations we present here consists in shooting a scalar 
soliton against a dipole-mode vector soliton. All the sim- 
ulations discussed here have been performed using a split- 
step operator technique using FFT, with grid sizes of 
up to 512x512 points covering a rectangular domain of 
68x34 adimensional units. The initial data are always a 
combination of stationary states. For instance, when a 
soliton is launched against a dipole, we start with 

u(x, 0) = u dipo i e (yL) + u so iiton(x - d)e~ lp ° x , (2a) 

t>(x,0) = V d ipole{x)- (2b) 

Here d = (d x ,d y ), d x 3> d v , d y is the impact parameter, 
and po is proportional to the initial (linear) momentum 
of the incoming scalar soliton. The initial data Udipoie, 
UsoMton, and v dipole are obtained numerically by a suit- 
able minimization procedure outlined in Ref. ||. 
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FIG. 1. Soliton-dipole scattering, (a) Snapshots of the in- 
tensity profile of each beam, (b) 3D plot of the total intensity 
|it| 2 + \v\' 2 , which shows the rotation induced in the dipole. 
(c) Same as (b), but with separated u and v components. 

The result is an inelastic collision in which the soli- 
ton becomes deflected and the dipole gains both linear 
and angular momenta. The whole process is depicted in 
Fig. [l]. Soliton scattering occurs when the incident beam 
has medium to large linear momentum or when it has an 
appropriate initial phase. For instance, in Fig. [l| the inci- 
dent soliton has sign (— ) and it crashes against the part 
of the dipole with (+) sign. A conservation law forces 
the dipole to rotate and the soliton becomes deflected, 
sometimes as much as by a 90 degrees angle. 
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FIG. 2. (a, b) Components of the linear momentum of 
the incident soliton (solid line) and dipole (marked by circles) 
after an inelastic collision with a large incident momentum, 
p u = J u'Vu dri, as a function of the impact parameter d y , 
which shows the crucial role of the dipole asymmetry. Total 
p y is not zero because of resultant radiation which is not seen 
in the figure. 
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FIG. 3. Absorption of a soliton by a dipole: (a) an inten- 
sity profile of each light beam -the darker the more intense; 
(b) 3D plot of the total beam intensity; (c) 3D plots where 
both u and v components have been separated. 

When the linear momentum of the incident soliton is 
large, it moves too fast to suffer a destructive influence 
from the dipole. In Fig. || we plot the exchange of the 
linear momentum between the soliton and dipole as a 
function of the impact parameter. The effective interac- 
tion is clearly attractive: the soliton coming from below 
(d y < 0) feels the drag of the dipole above it and gets 
deflected upwards (p y > 0), while the dipole moves down- 
wards. 

Soliton absortion by a dipole. The second family of nu- 
merical experiments is performed with solitons which are 
slow and, as is usual in scattering processes, the effects of 
the interaction process may be more drastical. For some 
impact parameters the soliton gets too close to the lobe 
of the dipole with the smallest phase difference and fuses 



with it with some emission of radiation and a subsequent 
rotation of the dipole. This is well reflected in Fig. || 
(radiation is not seen). 

Dipole-dipole collisions. The third family of numerical 
simulations corresponds to shooting dipoles against each 
other. These collisions provide a rich source of phenom- 
ena depending on the mutual orientation of the dipoles 
and on the initial energy. Figure |] summarizes the main 
results observed. There we see three cases (a-c) in which 
the dipole solitons are preserved. The figure shows an 
in-phase collision with weak interaction [Fig. 3(a)], an 
out-of-phase collision with repulsion [Fig. ||(b)], and an 
example of the collision with nonzero impact parameter 
in which two vortex states are created and they decay 
into a pair of spiralling solitons [Fig. |](c)]. 
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FIG. 4. 


Collisions of two dipoles 


against each other, seen 



from the center of mass of the system. 

The last case, Fig. ||(d) , shows an interesting inelastic 
process when two dipoles fuse into a more complex state 
which then decays creating a new dipole and a pair of 
simple solitons. All these processes may be understood in 
terms of the phase of the lobes of each dipole as described 
above. 

Experimental results. Generation of an isolated dipolc- 
mode vector soliton was reported earlier in Ref . Q . The 
dipole-mode soliton was created using two different pro- 
cesses: (i) phase imprinting, when one of the beam com- 
ponents is sent through a phase mask in order to imprint 
the required phase structure, and (ii) symmetry-breaking 
instability of a vortex-mode composite soliton. In this 
way, we obtain a dipole-like structure with a phase jump 
along its transverse direction that is perpendicular to the 
optical axis of the crystal [see Fig. ||(a), the beam u]. 
That dipole-like beam is then combined with the second, 
nodeless beam [the beam v in Fig. ||(a)], and the result- 
ing composite beam is focused into the input face of the 
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photorefractive SBN crystal (the crystal has the same 
parameters as in [Q), biased with the DC field of 1.5- 
2.5 kV applied along an optical axis. To control the de- 
gree of saturation, we illuminate the crystal with a wide 
beam derived from a while light source. Propagating 
in an effectively self-focusing saturable medium, such a 
composite input beam creates a dipole-mode vector soli- 
ton, as shown in Fig. |B|(b) (both components are shown 
separately). As discused above, the beam v creates an 
effective asymmetric waveguide that guides a dipole-like 
mode in the form of two out-of-phase solitary beams that 
mutually repell each other. 




FIG. 5. Experimental demonstration of the soliton-dipole 
collisions. (a,b) Formation of a dipole-mode soliton from two 
out-of-phase coherent beams (u) and co-propagating beam 
(v): (a) input beams, (b) output beams. Experimental pa- 
rameters are: V = 1.8 kV, z = 10 mm, and the initial powers 
are P v = 2 fiW, P u = 2.6 fiW. (c,d) Soliton-dipole interaction 
shown prior (c) and after (d) the collision, only the dipole 
component of the composite soliton is shown. Rotation of 
the dipole-mode vector soliton after the interaction is clearly 
visible. 

To observe the soliton-dipole interaction effects, we 
launch a scalar soliton beam against the dipole soliton. 
The input state is shown in Fig. |E|(c), where the dipole- 
mode soliton is presented by its two-lobe u component 
only. When the soliton interacts with a dipole, it gets 
deflected and transforms a part of its linear momentum 
into an angular momentum of the dipole that starts ro- 
tating is visible in Fig. |B|(d). A qualitative comparison 
between the theory and experiment is hard to carry out 
since the original theoretical model is isotropic, while 
the bias photorefractive crystal is known to possess an 
anisotropic nonlocal nonlinearity || making the dipole 
rotation anisotropic, since the dipole structure along and 
perpendicular to the optical axis is different. 

Conclusions. We have studied the phenomenology of 
collisions of the recently discovered dipole-mode vector 
solitons with other nonlinear localized structures. Apart 
from checking the robustness of the dipoles against strong 
interactions, we have shown that in many cases they be- 
have qualitatively as tightly bound molecules of light, 
with two major degrees of freedom (rotation as a whole 
and oscillation of the lobes of the dipole) which can be 



excited by collisions. Sometimes the dipole excitation 
is so strong that the structure behaves as a pair of spi- 
raling beams earlier analyzed in Ref. ||. This is only 
one of many interesting phenomena observed in simula- 
tions which also include excitation of rotational modes 
by collision with a scalar soliton, annihilation or strong 
deflection of the incident soliton, etc. Even richer phe- 
nomenology is observed when two dipole-mode solitons 
are made to collide. It is remarkable that the rich be- 
havior observed here may be understood qualitatively in 
terms of the structure of the objects colliding and the rel- 
ative phases of the dipole components. Finally, we have 
also verified experimentally some of our predictions. 

A rich variety of the effects described here might make 
these objects good candidates for practical applicability 
in the field of integrated optics. In this sense, the dipole- 
mode vector soliton resembles the operation of an elec- 
tronic transistor since it is an asymmetric object whose 
response is nonlinear and depends on the directionality 
of the input. We think the behavior described here, in 
addition to its fundamental interest, may have wide ap- 
plicability in the future. 
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